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We present a new model of “Stealth Dark Matter”: a composite baryonic scalar of an SU(Nd) 
strongly-coupled theory with even Nn > 4. All mass scales are technically natural, and dark 
matter stability is automatic without imposing an additional discrete or global symmetry. Con¬ 
stituent fermions transform in vector-like representations of the electroweak group that permit 
both electroweak-breaking and electroweak-preserving mass terms. This gives a tunable coupling of 
stealth dark matter to the Higgs boson independent of the dark matter mass itself. We specialize to 
517(4), and investigate the constraints on the model from dark meson decay, electroweak precision 
measurements, basic collider limits, and spin-independent direct detection scattering through Higgs 
exchange. We exploit our earlier lattice simulations that determined the composite spectrum as well 
as the effective Higgs coupling of stealth dark matter in order to place bounds from direct detection, 
excluding constituent fermions with dominantly electroweak-breaking masses. A lower bound on the 
dark baryon mass 300 GeV is obtained from the indirect requirement that the lightest dark 

meson not be observable at LEP II. We briefly survey some intriguing properties of stealth dark 
matter that are worthy of future study, including: collider studies of dark meson production and 
decay; indirect detection signals from annihilation; relic abundance estimates for both symmetric 
and asymmetric mechanisms; and direct detection through electromagnetic polarizability, a detailed 
study of which will appear in a companion paper. 

PACS numbers: 12.60.-i, 95.35.+d, 11.15.Ha 

I. INTRODUCTION 

Composite dark matter, made up of electroweak- 
charged constituents provides a straightforward mech¬ 
anism for obtaining viable electrically-neutral particle 
dark matter that can yield the correct cosmological abun¬ 
dance while surviving direct and indirect detection search 
limits, e.g., W- In this paradigm, the dark sector con¬ 
sists of fermions that transform under the electroweak 
group and a new, strongly-coupled non-Abelian dark 
force. This was considered long ago in the context of 
technicolor theories, where the strong dynamics was do¬ 
ing double-duty to both break electroweak symmetry and 
provide a dark matter candidate m- 

In this paper, electroweak symmetry breaking is ac¬ 
complished through the weakly-coupled Standard Model 
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Higgs mechanism, while the new strongly-coupled sector 
is reserved solely for providing a viable dark matter can¬ 
didate. This dark sector is not easy to detect in dark 
matter detection experiments or in collider experiments, 
and so we give it the name “Stealth Dark Matter”. Ear¬ 
lier work in this direction includes , and except for 

[25H36j . was often limited by the inability to perturba- 
tively calculate the spectrum and form factors due to 
strong coupling. 

The proposed dark matter candidate is a scalar baryon 
of SU(N d ), and hence N D must be evenQ We take the 
dark fermions to be in vector-like representations of the 
electroweak group. Hence, the constituent dark fermions 
can acquire bare mass terms (fermion masses that do 
not require electroweak symmetry breaking) while also 
permitting Yukawa interactions that marry dark fermion 


1 Fermionic baryons arising from odd Nn were considered in 

Ref. I3ll where the limit M -J 10 TeV was found to avoid the 

direct detection constraints from the magnetic dipole interaction. 




electroweak doublets with singlets. This yields a theory 
in which dark matter couples to the Higgs boson in a tun¬ 
able way that is essentially independent of the dark mat¬ 
ter mass itself. This is somewhat analogous to dark sec¬ 
tor models with a dark U(l) portal (e.g., ITT7UTTI ). where 
the coupling to the Standard Model is tunable through 
an otherwise arbitrary parameter - the kinetic mixing 
between the dark U(l) and hypercharge. 

The existence of both electroweak-preserving and 
electroweak-breaking masses for the dark fermions pro¬ 
vides two main benefits. First, given that the Higgs bo¬ 
son couples electroweak doublets with singlets, the global 
flavor symmetries of the dark fermions can be completely 
broken to just dark baryon number. All mesons can 
decay through an electroweak process (e.g., electrically 
charged mesons through W exchange) or through the 
usual chiral anomaly (e.g., the lightest neutral meson). 
Ensuring that these particles decay before big-bang nu¬ 
cleosynthesis sets a weak lower bound on the Higgs inter¬ 
action strength. (This is in contrast with [2UJ [42] where 
additional interactions were required to ensure mesons 
decay, e.g., through higher dimensional operators). The 
second reason is related to the orientation of the chiral 
condensate after the dark force confines. Large vector¬ 
like masses for the dark fermions ensure that the conden¬ 
sate can be aligned toward the electroweak-preserving 
direction, and thus the dark sector leads to only small 
corrections to electroweak precision measurements. We 
estimate the size of these corrections in this paper. 

There are many appealing features of an electroweak- 
neutral composite dark matter candidate made up from 
fermions transforming under the electroweak group, in¬ 
cluding: 

• All of the dimensionful scales are technically nat¬ 
ural, since they arise from fermion masses (vector- 
like and electroweak breaking) and the confinement 
of a strong-coupled dark force. 

• Dark matter stability is an automatic consequence 
of dark baryon number conservation. No addi¬ 
tional global discrete or continuous symmetries are 
required. For Njj > 3, operators involving dark 
baryon decay are necessarily dimension-6 or higher, 
and thus safe from GUT-scale or Planck-scale sup¬ 
pressed violations of dark baryon number. 

• There are no dimension-4 interactions of the com¬ 
posite dark matter particle with the Standard 
Model except with the Higgs boson. The direct 
detection scattering cross section is thus automat¬ 
ically suppressed compared with a generic elemen¬ 
tary WIMP candidate. 

• Higher dimensional interactions of the dark mat¬ 
ter with the Standard Model are suppressed, in the 
nonrelativistic limit, by several powers of the dark 
matter mass. For a composite scalar, the lead¬ 
ing operators are charge radius (at dimension-6) 
and polarizability (at dimension-7). The impact 


of these (and other) operators on the dark matter 
scattering cross section in direct detection experi¬ 
ments has been studied in pGniiaiiM. 

• Interactions of the dark baryon through the neu¬ 
tral weak current, the charge radius interaction, as 
well as the contributions to the electroweak pre¬ 
cision T parameter, are simultaneously eliminated 
if the fermion interactions obey a global custodial 
5/7(2) symmetry. Additionally, as we will see, dark 
matter electric neutrality also follows from custo¬ 
dial 5/7(2). To simplify our analysis, here we will 
primarily study the subset of stealth dark matter 
parameter space in which the custodial 5/7(2) is 
preserved. (This simplification is very familiar from 
composite Higgs theories, e.g. m)- 

• The abundance of a strongly-coupled dark scalar 
baryon could arise through several mechanisms: 
an asymmetric abundance (such as through elec¬ 
troweak sphalerons mi or other mechanisms 1521 ), 
when the mass is not too large few TeV, or a sym¬ 
metric abundance, when the mass is large (perhaps 
~ C>(100) TeV) pniESImj . 

We focus mainly on a confining 5/7(4) gauge theory 
dark sector with dark fermions transforming non-trivially 
under the electroweak group. We apply our recent re¬ 
sults (35! using lattice simulations for the spectrum and 
effective Higgs interaction for 5/7(4). As emphasized in 
EH 135], this theory is well suited for lattice calculations 
since we are not interested in the chiral limit of vanish¬ 
ing dark fermion masses. Indeed, lattice simulations can 
efficiently simulate the parameter region where the dark 
fermion masses are comparable to the confinement scale, 
exactly where the perturbative estimates are least useful. 

The organization of the paper is as follows. In Sec. |H] 
we discuss the assumptions and requirements to con¬ 
struct our stealth dark matter model. In Sec. lIIIl we detail 
the dark fermion interactions and masses. In addition, 
we write the electroweak currents in terms of the dark 
fermion mass eigenstates of the theory, detailed in Ap¬ 
pendix [A] Until this point, the discussion of the model is 
general. In Sec. |IV[ we simplify the parameter space for 
phenomenological and calculational purposes, applying 
a global custodial 5/7(2) symmetry and taking the ap¬ 
proximately symmetric dark fermion mass matrix limit. 
Then in Sec. [V]we discuss the light non-singlet mesons in 
the theory, in particular their decay rates and constraints 
from non-observation at LEP II. In Sec. IVII we discuss 
the stealth dark matter contributions to the 5 parameter, 
and demonstrate the parametric suppression that hap¬ 
pens in several regimes. In Sec. |VII| we obtain the Higgs 
boson coupling to the dark fermions. Then in Sec. |VIII| 
we apply our previous model-independent results on the 
5/7(4) spectrum and effective Higgs coupling to stealth 
dark matter. We obtain the bounds on the parameter 
space from the non-observation of a spin-independent di¬ 
rect detection signal at LUX. We briefly discuss the relic 
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abundance of stealth dark matter in Sec. |IX| Finally we 
conclude with a discussion in Sec. [Xj 

II. CONSTRUCTING A VIABLE MODEL 
A. Basic assumptions 

We assume that the dark matter candidate is a com¬ 
posite particle of a non-Abelian, confining gauge the¬ 
ory based on the group SU(Nd) with Nf flavors of 
fermions transforming in the fundamental representation. 
The number Nf is restricted by only the condition of 
confinement. For reasons outlined in the introduction 
(abundance, detectability), the dark fermions carry elec- 
troweak charges. Our model includes a tunable Higgs 
“portal” coupling between the dark sector and the Stan¬ 
dard Model via dimension-4 Higgs couplings^ We do not 
consider QCD-colored dark fermions since with JVjW 3, 
dark baryons would not generally be color singlets 

B. Requirements 

We require dark matter stability to be automatic, aris¬ 
ing from a global symmetry. This motivates consider¬ 
ing the dark baryon of the non-Abelian dark sector to 
be the dark matter HHS]. In the presence of GUT-scale 
or Planck-scale suppressed operators, the stability of the 
dark baryon should be sufficient to avoid cosmological 
constraints. 

The requirement of a sufficiently preserved acciden¬ 
tal baryon number disfavors a dark 517(2) group. First, 
there is no automatic baryon number in SU( 2 ) because 
there is no fundamental distinction between mesons and 
baryons. Imposing a global U( 1) baryon number is pos¬ 
sible (e.g. see CD) but in addition baryon number vio¬ 
lating dimension-5 Planck-suppressed operators such as 
/dark/dark H / M r] must be absent, where / dark is the 
dark fermion. (Otherwise, the dark 517(2) baryon would 
decay on a timescale much shorter than the age of the 
Universe.) 

For Njj > 3, operators involving dark baryon de¬ 
cay are necessarily dimension -6 or higher, and thus safe 
from GUT-scale or Planck-scale suppressed violations 
of dark baryon number. SU(Nd) with odd Ap is a 
perfectly interesting theory, having been studied be¬ 
fore for Ap = 3 by our collaboration [31]. There it 
was found that a fermionic dark baryon has a magnetic 


2 Other portals, such as a dark gauged U(1) group that kinetically 
mixes with hypercharge, are neither present nor required here. 

3 The obvious exception, when JVp = N c = 3, is discussed in 
l3ll l55l , which is not a focus for us due to the baryons being 

fermions. Construction of QCD-singlet dark baryons with Njy = 
6,12,18,... may be possible, but we do not study this possibility 
further here. 


dipole interaction that leads to a significant contribu¬ 
tion to spin-independent scattering. Constraints from 
the XENONIOO experiment were satisfied only when the 
dark matter mass M > 10 TeV m■ This strong con¬ 
straint on the mass scale implies the model is difficult to 
test at near-future colliders. 

The magnetic dipole interaction (and other higher di¬ 
mensional operators that require spin) are absent when 
the dark baryon is a scalar. We are thus naturally led 
to SU(Njj) with even Njj > 4, for which the otherwise 
strong constraints from direct detection are weakened, 
lowering the scales of interest into a regime that can be 
probed by colliders and other detection strategies. 

We assume the dark fermions have masses Mf on the 
order of the SU(Njj) confinement scale A£>. If the masses 
were much smaller, the dark sector would contain light 
pseudo-Goldstone pions that transform under the elec- 
troweak group, which are strongly constrained by col¬ 
lider experiments. A dark sector with purely vector¬ 
like fermion masses has approximately stable electrically- 
charged mesons due to dark flavor symmetries. Con¬ 
versely, a dark sector with purely electroweak breaking 
fermion masses has a dark matter candidate that is ruled 
out by spin-independent direct detection through single 
Higgs exchange. (For example, quirky dark matter [T 6 ] 
is now completely ruled out by Higgs exchange, given 
the direct detection bounds from LUX [3] combined with 
the relatively light Higgs mass [Ml [ST].) Fermions with 
both vector-like and (small) electroweak breaking contri¬ 
butions to their masses can avoid both problems. 

We require the lightest dark baryon to be electrically 
neutral. We also require Higgs couplings at dimension-4 
to pairs of dark fermions. These two requirements im¬ 
pose restrictions on the electroweak charges of the dark 
fermions. 

One solution is familiar from old technicolor theories 
(e.g. [Mi M]): requiring the dark fermion charges to 
roughly satisfy |F| |T 3 | where T 3 is the 517(2 )l isospin. 

Choosing doublets (IX 3 1 = 1/2 under 517( 2)l) then gives 
a finite number of discrete possibilities. 

A simple model that satisfies all of these requirements 
is shown in Table [T[ The electric charges of the dark 
fermions in the broken electroweak phase are Q = ±1/2, 
ensuring all hadrons have integer electric charges. So long 
as the lightest Q = 1/2 and Q = —1/2 dark fermions are 
close in mass, the lightest baryon will be a scalar and elec¬ 
trically neutral. Finally, with the assignments shown in 
Table |TJ all gauge (and global) anomalies vanish, which is 
automatic with fermions that transform under vector-like 
representations of the SU(Njy) and electroweak groups. 


III. DARK FERMION INTERACTIONS AND 
MASSES 

The fermions F“' d transform under a global 17(4) x 
17(4) flavor symmetry with [517(2) x 17(1)] 2 3 4 surviving af¬ 
ter the weak gauging of the electroweak symmetry. From 
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Field 

SU(N d ) 

(SU(2)l, Y 

/ F?\ 



*-U) 

N 

(2,0) 

'■■(J) 

N 

(2,0) 


N 

(1, +1/2) 

F$ 

N 

(1,-1/2) 

f 4 

N 

(1, +1/2) 

f£ 

N 

(1,-1/2) 


Q 



+ 1/2 
- 1/2 
+ 1/2 
- 1/2 


SU(2)l- The mass terms explicitly break the re¬ 

maining [SU( 2) x C7(l )] 2 down to 17(1) x 17(1) where one 
of the 17(l)’s is identified with U (l)y. (In the special case 
when Mg 4 = M 34 , the global symmetry is enhanced to 
517(2) x 17(1), where the global 517(2) acts as a custodial 
symmetry.) Thus, after weakly gauging the electroweak 
symmetry and writing arbitrary vector-like mass terms, 
the unbroken flavor symmetry is 17(1) x 17(1). 

Electroweak symmetry breaking mass terms arise from 
coupling to the Higgs field H that we take to be in the 
(2, +1/2) representation. They are given by 


TABLE I. Dark fermion particle content of the stealth dark 
matter model. All fields are two-component (Weyl) spinors. 
SU(2)l refers to the Standard Model electroweak gauge 
group, and Y is the hypercharge. In the broken phase of 
the electroweak theory, the dark fermions have the electric 
charge Q = T3 + Y as shown. 


this large global symmetry, one 517(2) (diagonal) sub¬ 
group will be identified with SU(2)l, one 17(1) subgroup 
will be identified with U(l)y, and one 17(1) will be iden¬ 
tified with dark baryon number. The total fermionic con¬ 
tent of the model is therefore 8 Weyl fermions that pair 
up to become 4 Dirac fermions in the fundamental or 
anti-fundamental representation of SU(Nd) with electric 
charges of Q = T 3yL + Y = ±1/2. We use the notation 
where the superscript u or d (as in F u , F d and later ip u , 
tp d , 4'“, \I/ d ) denotes a fermion with electric charge of 
<5 = 1/2 or Q = — 1/2 respectively. 

The fermion kinetic terms in the Lagrangian are given 
by 

£ D £ E il-f ""l■ ( 1 ) 

i= 1,2 i=3,4 m ,j=u,d 

where the covariant derivatives are 

Dy, = d„- igW£cr a /2 - ig D G\t h (2) 

D 2 =d^~ igW£a a /2 + ig D G b ^t b * (3) 

K, = d ,- - ig D G\\t b (4) 

D{ ll = d il -ig'Y^B lx + ig D G b tl t h * ( 5 ) 

with the interactions among the electroweak group and 
the new SU(Np). Here Y u = 1/2, Y d = —1/2 and t b 
are the representation matrices for the fundamental of 
SU(N d ). 

The vector-like mass terms allowed by the gauge sym¬ 
metries are 

C D M 12 e l0 F{Fl - M 3A F 3 F d + M d 4 F d F? + h.c., (6) 


C D y^eijFiWFt + y d u F 3 ■ rfF? 

- y&eijFi&F* - y£ 3 F 2 ■ H^F 3 + h.c ., (7) 

where again the relative minus signs are chosen for later 
convenience. After electroweak symmetry breaking, H = 
(0 u/\/ 2) t , with v — 246 GeV. Replacing the Higgs 
field by its VEV in Eq. 0 , we obtain mass terms for the 
fermions, in 2 -component notation, 


C Z> — (F“ F£)M U 
+ h.c ., 



(F d F d )M d 



( 8 ) 


with the mass matrices given by 


( M 12 y\ J 4 v/V2\ 

\V^v/V2 M 3i J 
_ ( M 12 yf 4 v/Vl\ 

\y23 v /^2 m 34 J 


(9) 

( 10 ) 


These Yukawa couplings break the remaining 1/(1) x 
1 /( 1 ) flavor symmetry to U(1)d dark baryon number. 
The mass matrices M u and M d correspond to the masses 
of two sets of fermions with electric charge Q = +1/2 
and Q = —1/2 respectively, in the fundamental repre¬ 
sentation of SU(Nd). The two biunitary mass matrices 
can be diagonalized by four independent rotation angles 

( M q m 2 ) = R wr lM * R ffi) (ii) 

(f M d ) = R ^y lMdR ( 9 ") - (I 2 ) 


where e 32 = e u d = — 1 = —e 12 and the relative minus 
signs between the mass terms have been chosen for later 
convenience. The mass term M\ 2 explicitly breaks an 
[517(2) x Z7(l )] 2 global symmetry down to the diagonal 
SU( 2)dia g x U( 1) where the 5f7(2)di ag is identified with 


where the rotation matrices are defined by 


R(o{) = 


( cos Of 
sin Of 



(13) 
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The 2-component mass eigenstate spinors are thus 


IV. SIMPLIFICATIONS 





t;) = m) 



= iR(0f) 


= iR(.o d 2 ) 



(14) 

(15) 

(16) 
(17) 


where the extra phase in Eqs. ensures the Q = 

— 1/2 fermions will have positive mass eigenvalues. 

The Lagrangian for the fermion mass eigenstates be¬ 
comes 

2 

£ A ~ E + Mftfx? + h.c.) (18) 

i— 1 

where the mass eigenvalues are M “ 2 for Q = 1/2, and the 
distinction between fermions if) and \ allows us to write 
the Q = —1/2 fermion masses as M-f 2 . The Dirac spinor 
mass eigenstates are constructed from the 2 -component 
Weyl spinor mass eigenstates in the usual way, 


T” = 


+/ = 


giving the Dirac fermion masses 


w \ 

V’fV 

i = 1,2 

(19) 

xt\ 

xf 

* = 1,2 

( 20 ) 


CD- 


A 

£( 

i =1 


—d. 




( 21 ) 


The fermion masses themselves are obtained from a 
straightforward diagonalization of the mass matrices, 

1/2 


Ml 2 = 


M\2 + IW 34 


-F 


m 12 - M 3 “ 4 




with mixing angles 
tan 20 “ = 


2V2v(M 12 y% 3 + M&yl 4 ) 


2M? 2 - 2(M 3 “ 4 )2 + (y“ 4 u)2 - (y 2 » 2 


( 22 ) 


(23) 


tan 20 n = 


2y / 2u(M 12 y/ 4 + M&yh) 

2M/ 2 - 2(M 34 ) 2 - (y“ 4 u) 2 + (y 23 ^) 2 


(24) 


with identical expressions for M 4 2 and tan 2 Of 2 with the 
replacement u o d everywhere. 

It is important to note that the electroweak currents 
(j+, jl, j 3 , jy) play an important role in the upcom¬ 
ing phenomenological discussions. Due to the extended 
expressions for these quantities in terms of our Dirac 
spinors, we have relegated a detailed derivation of the 
electroweak currents to Appendix [A] 


Our main interest is the more specialized case where 
the lightest Q = + 1/2 and Q = — 1/2 fermions are degen¬ 
erate in mass to a very good approximation. This leads 
to a neutral scalar baryon with a vanishing charge radius. 
While there are several ways this could be accomplished, 
we can simply impose a custodial 527(2) global symme¬ 
try on the Lagrangian. In order to simplify notation, we 
define c? = cos#?, sj = sin#/ and P L ^ R = (1 + 7 5 )/ 2 . 
In the custodial SU( 2) symmetric theory, c“ = cf and 


ST = sT 


A. Custodial SU(2) 


An exact custodial SU (2) symmetry implies the masses 
and interactions are symmetric with respect to the inter¬ 
change u ■O' d. This means the Lagrangian parameters 
satisfy 

2/14 = 2/14 — 2/14) 2/23 = 2/23 — 2/23) (25) 

m / 4 = m 34 = m 3 4 . 


Defining the overall vector-like mass scale M and differ¬ 
ence A to be0 


M = 


Mia + M 3 4 


A = 


M 12 — M 34 


(26) 


the dark fermion mass eigenvalues are 


My 2 =Ml|/i 2 + 2/142/23/^ _ (27) 

No u or d labels are necessary, since custodial SU(2) 
symmetry implies that there is one pair of Dirac fermions 
with electric charge Q = (+1/2, —1/2) with mass Mi (the 
lightest pair), as well as a second pair of Dirac fermions 
with electric charge Q = (+1/2,—1/2) with mass M 2 
(the heavier pair). The spectrum is illustrated in Fig. [l] 
In the limit 2 / 14 ) 2/23 "A 0, the fermions acquire 
purely vector-like masses, and thus the chiral conden¬ 
sate of the dark force is aligned to a purely electroweak- 
preserving direction. In order that the chiral conden¬ 
sate’s electroweak-preserving orientation is not signifi¬ 
cantly disrupted, we consider small electroweak breaking 
masses, 2/14^) 2/23^ <C M. 

This leaves two distinct regimes for the spectrum, de¬ 
pending on the relative sizes of i/yuy^v and A. 


B. Approximately symmetric mass matrices 

A second simplification, useful to analytically and nu¬ 
merically evaluate our results, is to take 2/14 — 2 / 23 - The 


4 We assume A < M , such that fermion masses remain positive, 
to avoid further fermion field rephasings. 
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V. LIGHT NON-SINGLET MESON 
PHENOMENOLOGY 


> k 


M^ d 

M 

M u,d 


t 

± v/A 2 + 2/i42/23^ 2 /2 

_ I _ 


Q = 


1 

2 



FIG. 1. Illustration of the fermion mass spectra considered 
in the paper. Four Dirac fermions (M>“, 'if, 'k(f, ^2) have 
masses (M“, Mf, , M $). The u ( d ) fermions have electric 
charge Q = +1/2 ( Q = —1/2); we assume an exact custodial 
SU( 2) global symmetry that ensures each Q = +1/2 fermion 
is accompanied by a Q = —1/2 fermion with equal mass as 
shown in the figure. If A < (A » yd/IdyiW) the 

mass splitting is dominated by electroweak breaking (preserv¬ 
ing) masses that we call the Linear (Quadratic) Case. See the 
text for details. 


mass matrices Eqs. TO > are approximately symmetric. 
Specifically, we can write 


yi4 = y + e v , 2/23 = y -e„, M < \y\■ (28) 

and expand in powers of e y . For example, the dark 
fermion masses become simply 


Mi )2 = M + 



(29) 


Theories with new fermions that transform under 
vector-like representations of the electroweak group 
generically have enlarged global flavor symmetries that 
can prevent decay of the lightest non-singlet mesons and 
baryons. In the case of dark baryons, this is a feature, 
providing the rationale for the stability of the lightest 
dark baryon of the theory. 

In the case of the lightest non-singlet mesons, this can 
be problematic, since some of these mesons carry electric 
charge)^] Stable integer charged mesons are strongly con¬ 
strained from collider searches as well as cosmology. One 
solution is to postulate additional higher dimensional op¬ 
erators that connect a dark fermion pair with a Standard 
Model fermion pair This must be carefully done 

to avoid also writing operators that violate the approx¬ 
imate global symmetries protecting the stability of the 
dark matter. In the stealth dark matter model, however, 
electroweak symmetry breaking can provide the source 
of global flavor symmetry breaking, leading to the de¬ 
cay of the lightest charged mesons. (We will not dis¬ 
cuss the lightest neutral mesons, but they are generically 
more difficult to produce in colliders, and they will decay 
through essentially the same mechanism as we describe 
for the charged mesons.) 

The lightest electrically charged mesons are composed 
dominantly of the dark fermion pairs II + = (H/ftl/") and 
n- = (tE'Y'E'i) ■ We can estimate the lightest meson life¬ 
time by generalizing pion decay of QCD to our model. 
The relevant matrix element is (see, e.g., [HTj ) 

<0|ilaxiall n± }=*/nlA (30) 


where /n is the “pion decay constant” associated with the 
dark force in this paper. The axial part of the electroweak 
current can be read off from the electroweak currents 


given in Eqs. (A5),(A6| 


axial ^ Caxial'l'^+S'E'? 


(31) 


where 


to leading order in 0(e y ). 

The distinct regimes are thus yv A and yv <C A. In 
the Linear Case j;r > A, electroweak symmetry break¬ 
ing is (dominantly) responsible for the mass splitting be¬ 
tween and 'S> 2 ' d - I n the Quadratic Case yv <C A, 
the splitting is dominantly attributed to the vector-like 
mass splitting A. As we shall see, the primary distinction 
between these two cases is in the Higgs coupling to the 
fermion mass eigenstates: proportional to y for the Lin¬ 
ear Case and y 2 for the Quadratic Case, hence the case 
names. A similar observation was also found in Ref. (SD] ■ 

From this point forward unless noted otherwise, we as¬ 
sume the fermion mass parameters satisfy an exact cus¬ 
todial SU(2) and the mass matrices are approximately 
symmetric. 


^axial — 


C 1 C l c 2 c 2 


72 


(32) 


and j't axial is identical upon u +>• d. In the custodial 
limit, Eq. (25), the axial coefficient is 

(2/14 - 2/I 3 )7 


^axial 


7 2 (8 M 2 + (y 14 - y 23 ) 2 v 2 )(8A 2 + (y 14 + y 2 3) 2 v 2 ) 

(33) 


5 We use the term “lightest mesons” and not “pions” since the 
would-be global symmetry that protects pion masses is com¬ 
pletely broken by the dark fermion vector-like masses. Never¬ 
theless, we use the symbol II to denote the corresponding fields. 
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Some insight can be gained using approximately symmet¬ 
ric mass matrices, Eq. (28). We then obtain 


Caxiai 


t y yv 2 

2M \/2A 2 + y 2 v 2 

e y v j 1 

2 M 1 yv/(\/2 A) 


(34) 


Linear Case 
Quadratic Case. 


The decay width can be obtained from pion decay of 
QCD by replacing V ud in the Standard Model with c ax i a i 
for the dark mesons. Since the charged dark mesons of 
this model are much heavier than the QCD pions, there 
are many possible decay modes. For a general decay into 
a Standard Model doublet (//'), assuming m/ m//, 
the decay width is 

/^2 / m 2 \ 

T(n+ -> //) = ^/nmjmncLiai ( 1 ^ ) ' ( ' 35 ' ) 


If run > mt + mb, the dominant decay mode is expected 
to be n+ —> tb, otherwise II + —>• t + v t and II + —> sc, 
with branching ratios of roughly 70% and 30% respec¬ 
tively. Note that the decay width has several enhance¬ 
ment factors relative to the QCD pion decay width 


r(n+ -> //) 

r(7T fl+V/j.) 


u axial 


\v ud \ 2 


/n 


tM (^) ( — ) (36) 


mn 


VI. CONTRIBUTIONS TO ELECTROWEAK 
PRECISION OBSERVABLES 


Stealth dark matter contains dark fermions that ac¬ 
quire electroweak symmetry breaking contributions to 
their masses. Consequently, there are contributions to 
the electroweak precision observables of the Standard 
Model, generally characterized by S and T (571 |BS]. In 
the custodial SU( 2) limit, Eq. (25), the contribution to 
T vanishes. There is a contribution to S, controllable 
through the relative size of the electroweak breaking and 
electroweak preserving masses of the dark fermions. 

The S parameter is defined in terms of momentum 
derivatives of current-current correlators EOT, 


S = 167tI1( 
d 
dq 2 


3Y'(0) 
167t 


(37) 


flV 


q^q v 


X^(q 2 ) 


q 2 = 0 


X^{q 2 ) = / d A x e 




(38) 


where the currents j£(x) and jy( x ) for the stealth dark 
matter model are defined in Eqs. (A7) and (A8). Af¬ 


ter some algebra and identifications of symmetric con¬ 
tractions, these definitions of the currents in terms of 
4-conrponent fermion fields lead to the current-current 
correlator. In the custodial limit, we obtain 


where for simplicity we have neglected kinematic suppres¬ 
sion. As an example, if /n — mn — v, we find the lightest 
charged dark mesons decay faster than QCD charged pi¬ 
ons so long as c ax i a i > 10 -8 . This is easy to satisfy with 
small Yukawa couplings and dark fermion masses at or 
beyond the electroweak scale. 

We can now make some comments about existing col¬ 
lider constraints on non-singlet mesons. The lightest 
charged mesons IU can be pair produced in particle 
colliders through the Drell-Yan process, and will decay 
through annihilation of the constituent fermions into a 
W boson. Because the Drell-Yan production is mediated 
by a photon and the mesons have unit electric charge, the 
production cross-section is substantial, leading to robust 
bounds from LEP-II. For charged states near the LEP-II 
energy threshold, the dominant decay mode is expected 
to be n+ —> t + v T as noted above. Reinterpreting the 
LEP-II bound from the pair production of supersymmet¬ 
ric partners to the tau (with the stau decaying into a tau 
and a nearly massless gravitino), we find mu > 86.6 GeV 
Stronger bounds from the LHC may be possible, 
although existing searches do not yet give any significant 
constraints on the charged mesons EDI; we briefly high¬ 
light the signals in the discussion. 

Using our lattice results from Ref. [33], we can trans¬ 
late the experimental bound on the mass of the pseu¬ 
doscalar meson into a bound on the baryon mass, 
m,B > 245,265,320 GeV when the ratio of the pseu¬ 
doscalar mass to the vector meson mass is mu/fn-y = 
0.77,0.70,0.55. 


2<#(*)#(0)> = ( U Gf L + 22 Gf L - 12 G^ - 21 Gf L ) 

, 2 2 (11 /-ifiv | 22 /-ifiv 12 21 s-ifii' \ 

-h C 2 S 2 V ^RR U RR U RR ^RRJ 

+ clsl ( n G- + 22 G- ) + c 2 s 2 rcz + 22 G- ) 

- c lC2SlS2 ( 12 Gf* + 12 G% l + 21 Gf* + 21 G% l ) , (39) 

where the connected contributions to the correlation 
functions are given by 


ij G^ B = {V?'fP A *ffi'fPB*?)\ 


connected 


(40) 


Here, A, B = L, R and the flavor indices i, j = 1, 2, where 
it is understood that the flavors labeled 2 have larger 
fermion masses than the flavors labeled 1. Since the u, d 
flavors have the same mass, the u and d labels are inter¬ 
changeable (i.e. everything is written in terms of the u 
flavors). 

We can obtain expressions for the mixing angle coeffi¬ 
cients. Like the case of light meson decay, if we consider 
an approximately symmetric mass matrix, with Yukawa 
couplings given by Eq. (28), all of the mixing angle coef¬ 


ficients are approximately equal to each other, differing 
only at first order in e y , i.e., 


22 22 22 22 
C 1 S 1 ~ C 2 S 2 ~ 0^2 — C 2 S 1 ~ CiC 2 SiS 2 


1 


2 2 
y v 


4 y 2 v 2 + 2A 2 


[1 + 0(e y )...] 


1 

- 4 x 


Linear Case 


y 2 v 2 /{ 2A 2 ) Quadratic Case. 


(41) 
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In the Linear Case, the mixing angles are approximately 
equal Ci ~ si ~ C 2 — S 2 — l/v2. In the Quadratic 
Case, all of the contributions to the S parameter are 
suppressed by {yv/ A) 2 . To calculate the S parameter in 
general requires lattice methods, paying close attention 
to the heavy-light splitting of the fermions, M 2 — Mi. To 
a first approximation we expect that in the limit of small 
mass splitting, M 2 — M\ -C M, 

U AB ~ U AB ~ U AB ~ U AB ~ U AB’ \^ Z ) 


This gives for the current-current correlator 


2 (33( x )Jy(°)) - [clsl + c^s 


2 

tv 
c y u 


_ 

2M 2 LRl 


2cic 2 sis 2 ] G^ r 


(43) 


where all of the Gll and Grr contributions self-cancel. 
Hence, we see that the contribution to the S parameter 
is suppressed asM>t) or e y -C 1, as expected. 


VII. FERMION COUPLINGS TO THE HIGGS 
BOSON 


In terms of the gauge-eigenstate fields, the interactions 
of the Higgs boson with the dark-sector fermions are, in 
matrix notation, 


h 

+ V^ 

+ h.c.. 


(Fi f£) 



(44) 


These matrices are not simultaneously diagonalizable 
with the mass matrices, Eqs. § ,(|To|). This means that 
the Higgs boson in general has off-diagonal, “dark flavor- 
changing” interactions with the mass eigenstate fields. 
Explicitly, we find in terms of the mixing angles 


C D 





„u u u I a u u u 
c l*2 i/14 "r *l c 2 i/23 
u u u u u u 

c l c 2 i/23 *1*2 i/14 


c^y^-s^yt 3 )(n) +[u ^ d) 

_ „u r u n.u _ r u Q u„.u / l v 7 

5 l c 2?/l4 C 1 s 2 2/23 J \^2 / 


(45) 


In the custodial 517(2) limit, we can drop the u and d labels since the Higgs coupling matrix is identical for both 
sets of fields. If we further take the limit of an approximately symmetric mass matrix, Eq. (281, the Higgs couplings 
simplify to 


C D 


yh 

M 2 — Mi 



vE- 2 


yv —\/2A 
— \/2A — yv 



(46) 


We observe both diagonal and off-diagonal Higgs cou¬ 
plings to the fermions. The off-diagonal dark flavor- 
changing interactions vanish in the limit A —► 0 and 
e y —> 0. In this limit an enhanced flavor symmetry among 
the fermions is restored, and the analogue of the GIM 
mechanism forbids such interactions at tree-level. The 
off-diagonal Higgs couplings lead to an inelastic scatter¬ 
ing cross section when a single Higgs is exchanged. This 
is highly suppressed unless the mass difference M 2 — Mi 
is near the (non-relativistic) kinetic energy of the dark 
matter in galaxy. Two off-diagonal Higgs couplings can 
be combined in a loop involving one heavier dark fermion 
and double Higgs exchange, but this is suppressed by the 
square of the Higgs couplings times a loop factor, as well 
as by the mass of the heavier fermions. 

The single Higgs coupling to the lightest fermions is 
finally 


(47) 


where 




y 2 v 

M 2 — M\ 


G(e y ) 


y 

V 2 

y 2 y 

2A 


Linear Case 


Quadratic Case. 


(48) 


(Note also that the single Higgs coupling to the heaviest 
fermions v k 2 is identical up to an overall sign.) Depend¬ 
ing on the relative size of yv and A, the Higgs boson cou¬ 
ples linearly or quadratically proportional to the Yukawa 
coupling y. The additional suppression of yv /A in the 
Quadratic Case will imply that spin independent scat¬ 
tering through single Higgs exchange can be significantly 
weaker when the mass difference between the lightest and 
heaviest fermions is dominated by the electroweak pre¬ 
serving mass A. 


C D ymh'b 1 T 1 











VIII. DIRECT DETECTION BOUNDS FROM 
HIGGS EXCHANGE 


In a previous paper [23], we determined the model- 
independent bounds on direct detection from Higgs ex¬ 
change for a scalar baryon of SU( 4). The model- 
independent result was expressed in terms of the effective 
Higgs coupling to the baryon 


m B A B ) (A fO 

9b =- af) . (49 

V J 

The first factor, the baryon mass m B (divided by the 
electroweak VEV), as well as the third factor 


AB) 

Jf 


m B 


Mi dm B 
m B dMi 


are extracted from our lattice results J5Hj|. The second 
factor 


v dM\(h) 
= Ml dh 


y/2 Ml 

(yy ) 2 

2Mi A 


Linear Case 
Quadratic Case 


(50) 

provides the effective coupling of the Higgs boson to the 
fermions (multiplied by v/M \), and we have evaluated 
the derivative for the two cases in our model. 

Unfortunately, we cannot directly apply our previous 
results on constraints in a-m B space to the parameters 
of the stealth dark matter model. This is because we 
do not know the dark fermion mass, Mi, independent of 
the lattice regularization scheme. We can, however, con¬ 
struct a regularization-independent parameter, the effec¬ 
tive Yukawa coupling y e ff, that is closely related to the 
model parameters: 


2/eff = 


m B 

V s/2Mi 


m B 

VsAMi 


Linear Case 
Quadratic Case. 


(51) 


The a parameter is therefore 


Linear Case 

(52) 

Quadratic Case. 

Recasting our previous constraints in a-m B space into 
y e fi-m B space, we can identify the region of parameter 
space that remains viable. The constraints for the Lin¬ 
ear Case are shown in Fig. [2] and the Quadratic Case 
in Fig. [3] In the top two plots for the respective fig¬ 
ures, the region above the LUX bounds represents the 
excluded parameter space for the model at a given dark 
matter mass ( m B ) and effective Yukawa coupling (j/ e ff)- 
The figures show a clear qualitative trend in how the 
predictions change as a function of dark matter mass. 


2/eff 

m B 

2 v2 
2/off T 

mi 


In particular, the cross-section is independent of ms for 
the Linear Case and inversely proportional to m B in the 
Quadratic Case. The bottom plots in Figs. |2|3| shows the 
maximum y e s allowed for a given dark matter mass. By 
increasing the splitting A between the vector-like mass 
terms, significantly more y e ff parameter space becomes 
available. 


IX. ABUNDANCE 

We now provide a brief discussion of the relic abun¬ 
dance of stealth dark matter. In the regime where the 
dark fermions have masses comparable to the confine¬ 
ment scale of the dark force, calculating the relic abun¬ 
dance is an intrinsically strongly-coupled calculation. 
Unfortunately, this calculational difficulty is not easily 
overcome with lattice simulations, due to the different 
initial and final states. Nevertheless, it is straightforward 
to see that the relic abundance can match the cosmolog¬ 
ical abundance through at least two distinct mechanisms 
that lead to two different mass scales for stealth dark 
matter. In this section we discuss obtaining the abun¬ 
dance of stealth dark matter through thermal freezeout, 
leading to a symmetric abundance of dark baryons and 
anti-baryons. Separately, we consider the possibility of 
an asymmetric abundance generated through electroweak 
sphalerons. 


A. Symmetric Abundance 

In the early universe at temperatures well above the 
confinement scale of the SXJ (4) dark gauge force, the dark 
fermions are in thermal equilibrium with the thermal 
bath through their electroweak interactions. As the uni¬ 
verse cools to temperatures below the confinement scale, 
the degrees of freedom change from dark fermions and 
gluons into the dark baryons and mesons of the low en¬ 
ergy description. Some of the dark mesons carry electric 
charge, and so the dark mesons remain in thermal equi¬ 
librium with the Standard Model quarks, leptons, and 
gauge fields. Since the dark baryons are strongly cou¬ 
pled to the dark mesons, they also are kept in thermal 
equilibrium. As the temperature of the universe falls well 
below the mass of the dark baryons, they annihilate into 
dark mesons that subsequently thermalize and decay (or 
decay then thermalize) into Standard Model particles. 
The symmetric abundance of dark baryons is therefore 
determined by the annihilation rate of dark baryons into 
dark mesons. 

The annihilation of dark baryons to dark mesons is 
a strongly coupled process. We expect B*B nn, 
B*B —>■ 3n, and B*B — )■ 4n, (and to possibly more 
mesons if kinematically allowed) to occur, but we do not 
know the dominant annihilation channel. If the 2-to-2 
process B*B — > nn dominates, one approach is to use 
partial wave unitarity to estimate the thermally averaged 
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FIG. 2. Constraints on the stealth dark matter model in 
the Linear Case of the model. The top and middle figures 
show the predicted values for the smallest and largest fermion 
mass explored in our simulations (corresponding to the pseu¬ 
doscalar to vector mass ratio mu/mv = 0.55,0.77) as well 
as LUX bounds. Various y e g values are plotted on the fig¬ 
ure, where y e g fts yms/Mi in this case. The dark grey re¬ 
gion is excluded by the LEP constraints on charged dark 
mesons. The bottom figure displays the maximum y e g al¬ 
lowed for a given dark matter mass. Each of the green curves 
represents a different fermion mass in the lattice calculation, 
mu/mv = 0.55, 0.7, 0.77 from top to bottom, and the bottom 
red curve is the result in the heavy fermion limit. 


annihilation rate [53] [55], 


(av) 


47r(v L ) 


( 53 ) 


where (v - 1 ) ~ 2.5 at freezeout [55]. Matching this cross 
section to the required thermal relic abundance yields 


FIG. 3. Same as Fig. [2] but for the Quadratic Case of the 
model. In this case, y B g ~ yms/s/MiA. 


tub ~ 100 TeV. An alternative approach is to use naive 
dimensional analysis [7DHZ21, which appears to lead to a 
larger dark matter mass. 

If the 2-to-3 or 2-to-4 processes dominate instead, the 
additional phase space and kinematic suppression low¬ 
ers the annihilation rate, and therefore lowers the scalar 
baryon mass needed to obtain the cosmological abun¬ 
dance. For recent work that has considered the ther¬ 
mal relic abundance in multibody processes, see [251 HI]. 
Suffice it to say a symmetric thermal abundance of dark 
baryons will match the cosmological abundance for a rela¬ 
tively large baryon mass that is of order tens to hundreds 
of TeV. 
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B. Asymmetric Abundance 


Early work on technibaryons demonstrated that 
strongly-coupled dark matter could arise from an asym¬ 
metric abundance BHS|. The main ingredient to obtain 
the correct cosmological abundance involved the elec- 
troweak sphaleron - the non-perturbative solution at 
finite temperature that allows for transitions between 
vacua with different^ B + L numbers)^] In the early uni¬ 
verse, at temperatures much larger than the electroweak 
scale, electroweak sphalerons are expected to violate one 
accidental global symmetry, B + L + D number, leaving 
B — L and B — D numbers unaffected mini. Here 
D number is proportional to the dark baryon number, 
with some appropriate normalization (for examples, see 

muni). 


Given a baryogenesis mechanism, the electroweak 
sphalerons redistribute baryon number into lepton num¬ 
ber and dark baryon number. As the universe cools, 
the mass of the technibaryon becomes larger than the 
temperature of the Universe. Eventually, the universe 
cools to the point where electroweak sphalerons “freeze 
out” and can no longer continue exchanging B 1 D 1 and 
L numbers. The residual abundance of dark baryons is 
p ~ rriBnB where the number density is proportional to 
exp[—ms/Tsph], where T sp h is the temperature at which 
sphaleron interactions shut off. 

If the baryon and dark baryon number densities are 
comparable, the would-be overabundance of dark matter 
(from niB O’ m nuc i e0 n) is compensated by the Boltzmann 
suppression. Very roughly, ms ~ 1-2 TeV is the natural 
mass scale that matches the cosmological abundance of 
dark matter |)J. A crucial component of the early tech¬ 
nibaryon papers 0H6] is that the technifermions were in 
a purely chiral representation of the electroweak group, 
like the fermions of the Standard Model. 

In stealth dark matter, given an early baryogenesis 
mechanism (or other analogous mechanism to generate 
an asymmetry in a globally conserved quantity Bami- 
180]). it is possible that electroweak sphalerons could also 
lead to the correct relic abundance of dark baryons con¬ 
sistent with cosmology. 

There is one critical difference from the early techni¬ 
color models (as well as the quirky dark matter model): 
The dark fermions in stealth dark matter have both 
vector-like and electroweak symmetry breaking masses. 
This leads to a suppression of the effectiveness of the elec¬ 
troweak sphalerons by a factor of a, c.f. Eq. (50), leading 
to a somewhat smaller stealth baryon mass to obtain 
the correct relic abundance compared with a technicolor 


6 In this section, B refers to baryon number and is to not be con- 
fused with the field defined earlier 

7 In addition, an asymmetric abundance could be generated 
through other mechanisms, see Ref. m, in which case the mass 
scales and parameters depend on the details of the particular 
mechanism. 


model (all other parameters equal). A more quantitative 
estimate is complicated by several factors: 

• Determining how the electroweak sphaleron redis¬ 
tributes the conserved global charges in the pres¬ 
ence of fermions that acquire both electroweak pre¬ 
serving and electroweak breaking masses. To the 
best of our knowledge, this calculation has never 
been done. 

• Determining the precise temperature at which elec¬ 
troweak sphalerons shut off, in the presence of both 
the Standard Model and stealth dark matter de¬ 
grees of freedom contributing to the thermal bath. 

• The baryogenesis mechanism itself, that determines 
the initial B — L and B — D numbers. 


Given the exponential suppression of the asymmetric 
abundance as the dark baryon mass is increased, it is 
clear that the upper bound on the dark baryon mass is 
nearly the same as the technibaryon calculation (updated 
to the current cosmological parameters), when stealth 
dark fermions have vector-like masses comparable to elec¬ 
troweak symmetry breaking masses. (This case is, how¬ 
ever, constrained by the S parameter, see Sec. VI). We 
can therefore anticipate that a range of stealth dark mat¬ 
ter masses will be viable, up to about a TeV. More precise 
predictions require further detailed investigation that is 
beyond the scope of this paper. 


X. DISCUSSION 

We have presented a concrete model, “stealth dark 
matter”, that is a composite baryonic scalar of a new 
SU(Nb) strongly-coupled confining gauge theory with 
dark fermions transforming under the electroweak group. 
Though the stealth dark matter model has a wide param¬ 
eter space, we focused on dark fermion masses that re¬ 
spect an exact custodial SU(2). Custodial SU( 2) implies 
the lightest bosonic baryonic composite is an electrically 
neutral scalar (and not a vector or spin-2) of the SU(No) 
dark spectrum, and in addition does not have a charge ra¬ 
dius. This yields an exceptionally “stealthy” dark matter 
candidate, with spin-independent direct detection scat¬ 
tering proceeding only through Higgs exchange (studied 
in this paper) and the polarizability interaction (stud¬ 
ied in our companion paper Ell). Custodial SU( 2) also 
allows for stealth dark matter to completely avoid the 
constraints from the T parameter. While contributions 
to the S parameter are present, they are suppressed by 
the ratio of the electroweak symmetry breaking mass- 
squared divided by a vector-like mass squared of the dark 
fermions. We also verified the lightest non-singlet mesons 
decay rapidly (so long as e y ^ 0), avoiding any cosmolog¬ 
ical issues with stable electrically-charged dark mesons. 

Specializing to the case of N D = 4, we then applied 
our earlier model-independent lattice results [333 to the 
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parameters of stealth dark matter, and obtained con¬ 
straints on the effective Higgs interaction. We find that 
the present LUX bound is able only to mildly constrain 
the Higgs coupling to stealth dark matter for relatively 
light dark baryons. Even weaker constraints arise when 
the effective Higgs interaction is quadratic in the Yukawa 
coupling, which is a natural possibility when the two 
pairs of dark fermions are split dominantly by vector-like 
masses, i.e., yv <C A. 

While we have considered many aspects of stealth dark 
matter, several avenues warrant further investigation: 

• Chiral symmetry forbids additive renormalization 
of the fermion masses; we have focused on the 
regime where the constituent fermion mass is com¬ 
parable to the confinement scale Mf ~ A/j, 
since this is best-suited for lattice simulations, ex¬ 
actly where analytic estimates are least useful. It 
would be interesting to consider a broader range of 
fermion masses relative the confinement scale, to 
understand the relative scaling of the Higgs inter¬ 
actions. 

• A more precise calculation of the S parameter is 
possible using lattice simulations for the relevant 
correlators. This would allow us to place numer¬ 
ical bounds on the parameters of the theory, that 
could be stronger than the bounds from the non¬ 
observation through direct detection. 

• We would like to unpack y e g [c.f. Eq. ©] and 
obtain constraints on the Yukawa couplings of the 
model. However, this requires translating the 
fermion masses from the lattice regularization into 
a continuum regularization. 

• Dark meson production and decay at the LHC is 
ripe for exploration. Dark meson pair production 
would proceed through off-shell EW gauge bosons, 
qq —>■ n+n-, qq —> n°n°, and qq' —f n ± n°. These 
could have spectacular signals at the LHC. Neutral 
mesons decay into fermion pairs and dibosons (ex¬ 
plored in other related models in J5DJ @31 ES2i [S3])- 
For charged dark mesons, with masses in the range 
m n ± ~ 90—180 GeV, the decay n + —> t + v t domi¬ 
nates, while for masses above this, n + — > tb is dom¬ 
inant. Charged pion pair production could there¬ 
fore lead to tbbt signals with the tb and bt pairs 
reconstructing to the same mass. To the best of 
our knowledge, this type of resonance search is not 
being performed at the LHC. 

• More insight into the thermal abundance of stealth 
dark matter, perhaps using lattice simulations, 
would help narrow the interesting mass range that 
matches cosmological data. 

• Asymmetric production of stealth dark matter 
seems very promising, but has several calculational 
obstacles to overcome to arrive at a quantitative 


relationship between the abundance and the other 
parameters of the theory. 

• If stealth dark matter has an asymmetric abun¬ 
dance, there are potential limits from neutron star 
lifetimes [85H5S] though the precise bounds depend 
sensitively on the equation of state of the neutron 
stars. 

• There are tantalizing signals of a 7 -ray excess be¬ 
tween about 1-10 GeV in the galactic center (see 
for example f87H9T] ). A recent analysis [92] sug¬ 
gests that this could arise from dark matter up to 
300 GeV. It is intriguing to consider the 7 -ray sig¬ 
nal spectrum that could arise from a symmetric 
abundance of stealth dark matter with annihila¬ 
tion into a multibody final state [!J3] with mixtures 
of four or more heavy fermions and multi-gauge 
bosons (from BB* ->nn...-> SM states). 

Finally there are broader model-building questions to 
consider. One is the choice of scales Mf ~ Ap that has 
been the focus of this work. This could arise dynamically. 
For example, if there are sufficient flavors in the SU(No) 
gauge theory such that it is approximately conformal at 
high energies, then as the theory is run down through the 
dark fermion mass scale Mf, the dark fermions integrate 
out, and confinement sets in at A^ ~ Mf. This is well 
known to occur for supersymmetric SU ( N ) theories in 
the conformal window that flow to confining theories once 
the number of flavors drops below Nf < 31V/2 [S5]. The 
origin of the vector-like masses of the fermions is also 
an interesting model-building puzzle. However, just as 
SM fermion masses are vector-like below the electroweak 
breaking scale, we can imagine dark fermion vector-like 
masses could be revealed as arising from dynamics that 
breaks the flavor symmetries of our dark fermions at some 
higher scale. 
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Appendix A: Weak Currents 


We examine the dark fermion contributions to the elec- 
troweak currents. In the gauge eigenstate basis, the cur¬ 
rents are 


£ = --^= (f^v^F? + F^a^F^j 

£ = - -L + F 2 d VF 2 “) 

ft = E ~ Ff^Ft) 

Z i=l,2 

ft = E ( F?'<TF? - . 


(Al) 

(A2) 

(A3) 

(A4) 


In the mass eigenstate basis given by Eqs. (141-( 17), the 
currents can be rewritten in terms of the 4-component 
Dirac fermions defined by Eqs. (19),(20). After some 
algebra, one obtains 


■u 

ft = 


■u 

3- = 


■U 

3 3 = 


a/2 


nr ( c u Ap l + c u 2 4 pr) n + nr ( s u lS ip L + s u 2 s d 2 p R ) n 


+nr ( rsip L + 4 s d 2 p R ) n + nr (44 pl+s%4p r ) n 


l r —t 


\/2 L 


nr ( 4c\p l + c d c u 2 p R ) n+nr (*?«? p L + s d 2 s u 2 p R ) n 


l r^— 


2 L 


■u 

Jy = 


+nr (c d nPL+44 Pr) n+nr (44 p l + s d 2 c u 2 p R ) n 
nr ((r) 2 p L + (4) 2 p r ) n + nr {{si?p L + {s u 2 fp R ) n 

- nr (( 4) 2 pl + (. 4?pr) - nr (^/pl + (s d 2 ) 2 p R ) n 
+ nr (c“s“Pl+ 4 s%p r ) n ± nr w$p L +8 $<%p r ) n 

- nr ( 44p l + 44pr) ^2 - nr (44 p l + 44pr) n 

nr ((4 ) 2 pl + {s 2 ) 2 pr) n + nr {nfp L + o 4 fp R ) n 

- nr (( s d fp L + ( s d 2 ) 2 p R ) n - nr (( 4?p L + ( 4?p r ) ^ 

- nr {c\ s \P L + c u 2 s u 2 p r ) n - nr n 4 P L +s^p R ) n 
+ nr ( 44p l + 44pr) ^2 + nr (44 p l + s d 4 p R ) n 


l rA- 


2 L 


(A5) 


(A6) 


(A7) 


(AS) 


where 4 = cos0-, s{ = sind/ and Pl.r = (1 =F 7s)/2 
are the left- and right-handed projectors. In general, the 
dark fermions contribute to both the vector and axial 
currents with strengths given by the mixing angles. It is 
easy to verify that the electromagnetic current, 


,-a* = F + 

Jem j 3 ' JY 

= {crnm±Qrtrn 


(A9) 


with Q U 'd = ±1/2, is consistent with a pure vector cou¬ 
pling of the dark fermions to the photon independent of 
mass mixing angles. 

Interestingly, if the mass matrices Eqs. (|9|),(fTo|) are 
symmetric, i.e., y \ 4 = y% 3 and yf 4 = y$ 3 , then just two 
mixing angles are required, i.e., d 4 = d 2 and 9 d = 9 d . In 
this case, the mixing angles factor out of the left-right 
gamma matrix structure, leaving all of the electroweak 
currents to be purely vector (with vanishing axial cur¬ 
rent). This is unlike the Standard Model, where the 
SU{2)l currents are purely V — A. The difference be- 
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tween this model and the Standard Model is the struc¬ 
ture of the dark fermion mass matrices that include both 
vector-like and electroweak symmetry breaking masses. 
It is also interesting to calculate the neutral current 

3z = is ~ sin2 e wj & n • (A10) 

For the neutral baryon state, 

(B\f z \B) ~ (All) 

+ J (K) 2 + (c£) 2 - (c d f - (c d f) 

+ J (-(ci) 2 + (4) 2 + ( cf ) 2 - (c d 2 f) . 

In the limit of zero momentum exchange ( Q 2 = 0), the 
vector form factor (B|'I'i 7 a ‘^'i \B) evaluates to 1, while 


(c ?) 2 + (c^-(c?) 2 -(c^ 


r d \ 2 


Custodial symmetry violation is therefore restricted by 
requiring the coupling of the lightest neutral baryon to 
the Z boson be small enough to have evaded direct de¬ 
tection. There are several limits in which this can occur: 


the axial-vector form factor ( 2 ?| 1 ® r i 7 M 7 5 \f , i|.B) for a scalar 
baryon vanishes. In the presence of an exact custodial 
SU(2) symmetry, which is the focus of this paper, we 
have cf = c d and the Z coupling vanishes identically at 
any momentum exchange. 


On the other hand, if custodial symmetry is broken, 
then the lightest neutral baryon acquires tree-level cou¬ 
plings to the Z boson. To illustrate the size of these 
couplings, consider taking the dark fermion mass matri¬ 
ces to be exactly symmetric (?/23 = Via) but allowing for 
a small, custodial symmetry-violating difference in the 
Yukawas, y u = y + £ and yd = y — g where £/y <C 1. The 
coefficient of the weak neutral vector current becomes 


2 \/ 2 - — 
y yv 

g ( w) 2 

y A 2 


Linear Case 
Quadratic Case. 


(A12) 


£/y <C 1 (any scenario), A /{yv) -C 1 (Linear Case), or 
(yv)/ A <C 1 (Quadratic Case). This suggests that mod¬ 
est custodial symmetry violation is possible but rather 
constrained. 
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